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ABSTRACT 


A method independent of physical modeling assumptions is presented to 
analyze high-speed flame photography and cylinder pressure measurements from a 
transparent-piston spark-ignition research engine. The method involves defin- 
ing characteristic quantities of the phenomena of flame propagation and com- 
bustion, and estimating their values from the experimental information. Using 
only the pressure information, the mass fraction curves are examined. A new 
empirical burning law is presented which well simulates such curves. Statis- 
tical data for the characteristic delay and burning angles show that cycle-to- 
cycle fractional variations are of the same order of magnitude for both angles 
(about 20 percent). Using only the photographic information on flame front 
contours, statistics are obtained for the apparent true ignition delay time. 

Enf lamed volume, area of the flame front and an average normal flame front 
speed are estimated as a function of time. From the combined analysis of 
cylinder pressure history and flame photography the entrainment speed u g 
(often called turbulent flame speed) is estimated and found to increase rapidly 
during the initial period of flame propagation. At a later time, when the mass 
fraction burnt is greater than about five percent, u g remains approximately 
constant at a value of the order of 10 m/s for the present operating conditions 
with cyclic fractional variations of about 20 percent. Comparison of the- en- 
flamed and burnt mass fractions indicates that a substantial amount of unburnt 
mass is inside the enflamed region during flame propagation. Comparison of the 
rates of entrainment and burning shows that the two processes are not coupled 
as it is implied when a thin flame is assumed. The characteristic burning time 
is found to be of the order of 1 ms and tends to decrease through the combustion 
period. The flame stretch factor is also estimated. The present work concen- 
trates on experimental evidence which can be used to test theoretical models 
and computer simulations of the phenomenon. 



1 


1. INTRODUCTION 

It is only in the last few years that the historically established empir- 
ical method of developing an internal combustion engine is being intensively 
and progressively complemented by the science of modelling. The reasons for 
this change are found in two main problems of the present society: the envir- 

onmental problem, with the new constraints of strict emission standards, and 
the fuel shortage problem, with the consequent rising costs of the primary 
energy sources and the demand of more efficient use. 

The empirical method of development cannot face, in an economic way, the 
complexity of the problem of optimizing efficiency at acceptable power levels 
subject to the new, rigid environmental and energy constraints. The present 
goal of the science of modelling is to provide a better analytical basis for 
such engine optimization effort. But there is also a less pragmatic reason 
that has received more and more attention in the history of the internal com- 
bustion engine: the need to understand the governing physical mechanisms of 

the processes occurring in a real operating engine. 

The latter task is obviously of a great complication, since it involves 
going beyond the state of the art in all thermal sciences: non-equilibrium 

thermodyramics and chemical kinetics, unsteady fluid-dynamics and turbulence- 
all coupled in one problem. 

A number of different approaches are being tried to tackle the modelling 
problem, and their success clearly depends on how the calculated results com- 
pare with experiments. The aim of the present work has been to extract experi- 
mental information on that small but important portion of the whole operating 
cycle of the spark-ignition engine where the combustion process takes place. 
This has been done by developing and applying a simple method to analyze the 
data generated by a non-conventional experiment--the non-conventional part 
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consisting in high-speed photography of 'flame propagation through a transparent 
piston in a single-cylinder, spark-ignition research engine. 

It is well known that turbulence plays a fundamental role in the opera- 
tion of a real engine. Mathematical theories have been proposed to deal .with 
the turbulence problem [1 - 6], but no established physical description of a 
turbulent flow is yet available. The physics of turbulent flame propagation 
is even less understood [7-19], For such reasons an attempt is made in this 
work -to define quantities and obtain experimental estimates in a way which 
is independent of physical models of interpretation. The aim of this- effort 
is to obtain experimental evidence against which the assumptions of any model 
of the phenomenon may be tested. 

After a brief description of the experiment, a model independent 'method 
of analysis of the experimental data is presented and the results- discussed. 
Consistent with the declared effort .of the present work, physical interpreta- 
tion of the results is left for future work. 

2. EXPERIMENTAL APPARATUS AND MEASUREMENTS 

The experiment was carried out on .a transparent-piston engine [44]. Figure 1 
shows an axial section of the cylinder-piston system which allows optical 
access to the combustion chamber through the quartz piston window. The cylin- 
der head is a section of the head used by Ford on their 400 CID V-8 engine. 

Table 1 gives the details of the engine geometry and valve timings 

Because optical access is obtained through the piston, the original cham- 
ber and valve configuration could be conserved'. The purpose of this experi- 
ment was to obtain detailed photographic information on the flame profiles 
through each cycle. High-speed photography was taken with a Hycam Model 
K-2001-R- rotating-prism camera capable'of taking up to five thousand frames 
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per second. In order to increase flame luminosity, the mixture of iso-octane 
and air was seeded with uniformly suspended salts, appropriately chosen not to 
affect flame propagation. The edge of the film was marked by a neon light 
pulse, synchronized to the flywheel, every twenty crank-angle degrees. Such 
markers provide a means to associate a crank angle to each frame of the photo- 
graphic sequence. 

The more conventional aspect of the experiment is the measurement of in- 
take mass flow rates and of the cylinder pressure history. For the pressure 
measurement, the output of a piezoelectric transducer, suitably amplified, was 
A/D-converted by an on-line POP 11 digital computer storing one pressure value 
every crank-angle degree. 

The experimental procedure was the following: first, the engine was 

motored for a sufficiently long time to reach a complete steady state. Then, 
an electronic switching system was activated for the following sequence of 
operations: start the high-speed camera*, wait for its speed to reach a thresh- 

old value; light the spark for ten to thirteen cycles, simultaneously activa- 
ting the pressure data acquisition system. 

As a result of this procedure the complete data set for each cycle of a 
sequence consists of the pressure measurement and the synchronized flame 
photographs as a function of crank-angle. Table 2 summarizes the operating 
conditions of the six sequences for which a complete data set was available. 

Figure 2 shows a typical sequence of pressure traces as a function of 
crank angle for a sequence of 10 successive firing cycles. A consequence of 
the experimental procedure is that the first cycle of each sequence has no 
residual burnt fraction x r and its pressure trace departs from the average. 

Figure 3 shows a typical photographic sequence corresponding to pressure 
trace number 9 in Figure 2. The outer circular crown of the piston is used 
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to retain the quartz window (see Figure 1) and no direct photographic informa- 
tion on flame profiles is obtained in the corresponding corner region of the 
combustion chamber. 

3. ANALYSIS OF EXPERIMENTAL DATA 
3.1. Mass Fraction Burnt 

Given a pressure measurement as a function of crank angle, a measure of 
inlet mass flow rates and an estimate of the residual burnt fraction x^, the 
energy balance equation for the combustion chamber can be solved to yield the 
mass fraction burnt = m^/m and all the other relevant thermodynamic proper- 
ties of the gas mixture at each instant of time through the cycle. To perform 
this task, a version of the thermodynamic model presented by Tabaczynski, 
Ferguson and Radhakrishnan [20] was used. The estimate of heat losses was made 
following the correlation of Woshni [21] and crankcase blow- by was neglected. 

From a different analytic version of the energy balance (Beretta [22]), it 
can be shown that neglecting blow-by will result in an underestimate (of the 

order of the blow-by itself) of the maximum value of x, , x. . Given a pres- 

b b max ^ 

sure trace, many factors may influence the estimate of x b and its maximum 

value, and they are all coupled in the energy balance equation. Underestimating 

■heat losses, mass blow-by and pressure, results in underestimating x^ and 

x b max' Underestimating inlet mass flow rates (i.e., the initial charge mass) 

results in overestimating x. and x. 

b b max 

Inlet mass flow rates were determined with a reasonable accuracy, and 
heat losses have a slight influence on x^ max > The problem of the dynamic 
calibration of the pressure transducer is, on the contrary, an important one. 
When the statically calibrated pressure trace was used to estimate x^, low 
values of x^ max were obtained (in the range 80 to 85 percent). To find what 
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reasonable value x. , should assume in our engine, the mass fraction quenched 
b max 

at the chamber walls was estimated to be about 1 percent of the charge mass 

(data from Ferguson and Keck [23] and Weiss [24] were used), the mass trapped in 

the piston crevice above the upper ring was about 4 percent, and a reasonable 

estimate of the blow-by underestimate was about 2 percent. The conclusion was 

that x. should be around 93 percent. The disagreement was attributed to 
b max r 

the calibration procedure, and it was decided to scale the pressure traces to 
obtain, on the average, the expected x b max - 

Obviously, having done this, a check had to be made on what is its 
influence on the subsequent analysis. When the results of the complete analy- 
sis for the same run, with and without a scaling factor of 10 percent were 
compared, the important quantities that will constitute the final results did 
not vary by the same amount as x, a ; namely, estimates of entrainment speed 
varied by a factor of 2 percent, well below the expected accuracy from such an 
experiment. Estimates of mass entrainment rate varied by less than 1 percent, 
while mass burning rates, as expected, varied by about the same 10 percent. 
Nonetheless, all trends and phase effects were conserved and we reached the 
conclusion that the scaling factor does not significantly affect the final 
results. 

As a result of the thermodynamic analysis the following quantities of 
importance to the subsequent analysis are known as a function of crank-angle: 

p, V, m pressure, volume and total charge mass 

p , p , unburnt density, its time-derivative and unburnt 

gas temperature 

x b , x , x b mass fraction burnt and its first and second time- 


derivatives 
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Figure 4 shows the normalized mass fraction burnt n = x,/x, vs. 0 

b b max 

for the pressure curves in Figure 2. 

In Figure 5, among other definitions, that of "burning angle" is visual- 
ized: 

A8 b - (1) 

x b max 

where N is the crankshaft speed (rev/min) and A6 b represents the number of 
crank-angle degrees that a line, tangent to the mass-fraction-burnt curve at 
the point of maximum slope, takes to rise from x b = 0 to x fa = x b max - The 
interval between spark and the angle at which the same line intersects - 0 
is called "delay angle", A0^. Table 3 shows the statistics of Ae^, AS^ and 
A0J + A0^ for the analyzed sequences. Contrary to the common interpretation 
that the variations in delay angles are the main cause of cyclic variations, 
it is observed that the relative variations of both delay and burning angles 
are of the same order (about 20 percent). That positive variations of AO^ 
correlate to -negative variations of AS^ can be inferred from the lower cyclic 
variations of the sum of the two angles. 

3.1.1. Empirical burning law 

In spite of the large cyclic variations in the mass-fraction-burnt curves, 

a remarkable stability is observed for the normalized plot of x,/x, vs. 

b b max 

x b^ x b max P resen ^ ec * in Figure 6 where the experimental curves for about seventy 
cycles have been superimposed. 

It is found that the following family of curves well represents such 
behavior: 

-1 

f(n) = C(a,b) [brf a + (1 - b) (1 - n)” a ] (2) 


where 
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n = x,/x, = normalized mass fraction burnt 

b b max 

f <"> = V*b max ■ " x b max /x b max = ^ 8 b /6N 
a = contact parameter 


b = symmetry parameter 


C(a,b) = b 
n(f = 1) = 



= normalization factor 


= n at the point of 
maximum slope 


Shown in the upper part of Figure 6 are plots of f(n) vs n for a = 0.84 and 
b = 0.45 - 0.10, the value of the parameters best fitting the data of the 
lower part of the figure. Equation (2) is easily integrated to yield 


e(n) - e s 



(3) 


where 6 s is the crank-angle at spark and g(n) is a dimensionless crank-angle 
(or time) 

g(n) = bp 1 a - (1 - b) (1 - n)"*“ a + 1 - b (4) 

Figure 7 shows the effect of parameters on the empirical burning law n vs g(n). 

The high degree of contact, presented by g(o) after spark, is a charac- 
teristic that allows the use of 9 s as the spark crank angle. In the other 
existing empirical burning laws (Wiebe [25]; Blumberg and Kummer [26]), 6 $ is 
not identified with the spark-angle, but as the point at which combustion 
begins, thus assigning to the ignition delay time a separate role, decoupled 
from flame propagation. This commonly accepted decoupling procedure is now 
believed not to represent the physics of the phenomenon. For example. Hires, 
Tabaczynski and Novak [27] have recognized, in their model, that the same, 
mechanism should govern both delay and burning periods of combustion. Only, 
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due to the nature of turbulence and to variations in the mean flow field 
around the spark plug at spark time, the initial flame kernel may find a dif- 
ferent geometrical and mean flow environment in which to develop. 

Equation (4) makes it possible to link the values of all the intervals 
that can be defined on the mass-fraction-burnt curve, to the value of A0 b pre- 
viously defined. With reference to Figure 5, the definitions of the other 


intervals are: 


> 

CD 

O 

• 

= 0(0.1) - 0 S 

A6 d 

= 0(n(f = 1)) - 6 S - j A0 b = delay interval 

ie 0.1-0.9 

= 0(0.9) - 0(0.1) 

> 

CD 

O 

= 0(0.99) - 0(0.01 ) = combustion interval 

A0 t 

0(0.99) _ 0 S ~ total combustion interval. 


In Table 4 the values of Ae^, A0 b , A0 Q -j , A0 Q g and A0^ are shown for the 
curves of Figure 4. In round brackets are their ratio to AO^. The ratio 
A0 t /A0 b , calculated from the burning law, is found to be very sensitive to the 
parameter a. For this reason it is kept fixed to the value 0.84 while the 
parameter b shown in Table 4 is calculated to match the experimental values of 
A0 d /A0 b . In square brackets are the corresponding predicted values for 
a = 0.84 and b = 0.42. For the given engine geometry, operating conditions 
and the case of edge ignition, the empirical burning law could be used to 
simulate cyclic variations by keeping the parameter a fixed and chosing the 
values of b and A0 b from Gaussians of given average and standard deviation (in 
our case both of the order of 20 percent). 

3.2. Flame Front Geometry 


In this work the flame front is operationally defined as the edge of the 
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photographic image on the film. Spatial averaging over small scale structure 
is already included due to the finite resolution of the system and line of 
sight integration of the exposing radiation. For each cycle of a sequence, 
the flame front profiles corresponding to each frame were superimposed on a 
drawing. The coordinates of the significative points for each flame contour 
were introduced in a digital computer for further automatic analysis as is 
shown in Figure 8a. 

To obtain complete geometrical information on flame contours in the 
whole combustion chamber, assumptions are required to reconstruct the non- 
visible regions. The simplest reasonable assumption, generally consistent 
with observations [28], is to fit to the digitalized points of each flame 
front profile the best (least squares) circle, leaving as free parameters its 
radius r^ and the position x , y c of its center. Thjs assumption is clearly 
less restrictive than holding the center fixed at the spark location. Figure 
8b shows the result of this procedure. The third dimension (i.e., the coord- 
inate z c of the best sphere) is assigned by assuming that the apparent flame 
center keeps the same relative distance from the cylinder head as the spark 
plug at spark time. Due to the high aspect ratio of the combustion chamber 
close to TDC, this assumption is not important. 

A consistency check of our assumptions is obtained by following the path 
of the best center with time. During the early period of flame propagation, 
the center is found to move out from the spark plug, along a circular path on 
the x, y plane, at a velocity of about 1 m/s, consistently with the swirl 
velocity observed on the films at the end of combustion. At later times, 
the best center normally moves outward from such circular path, consistently 
with the fact that the flame front profiles become flatter. 
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Having r f , x c , y , z Q and knowing the geometry of the combustion chamber 
at the corresponding crank-angle, it is a geometrical problem to intersect the 
best sphere and the chamber geometry and obtain the shape of the enflamed 
region. The result is the estimate of the following parameters of the enflamed 
region, which are important in the subsequent analysis: 

volume of the enflamed region 

V^-, volume of the non-enflamed region (V^, + = V) and its time 

derivative 

area of the flame front 

Ap^ area of piston surface not yet interested by the flame 
average normal speed of the flame front. 

A direct estimate of is obtained dividing the average distance of two 
flame iprofiles by the time between the corresponding frames. An integral 
estimate .of u^ is obtained from the knowledge of the above geometrical para- 
meters of the enfl aired region and the equation 

^f = "'“/f “ u p A pf ^ 

which defines the average normal speed of the flame front. 

Figure 9a shows the best flame front distance from the spark plug 
Tf + x c - x g versus crank-angle 0 for the typical cycle of Figure 8. Extra- 
polation at zero distance yields an apparent ignition crank-angle 6... Figure 
9b shows direct and integral estimates of u^; the remaining information in 
Figure 9, is discussed later. 
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3.3. Entrainment and Gas Speeds 

Figure 10a shows cylinder head geometry and assumed spherical flame fronts 
for an early and a late visible contours. Figure 10b shows a plot of time 
versus flame front distance from spark. The average particle path is also 
indicated. 

Assuming uniform pressure throughout the combustion chamber and negligible 
volume of the chemically reacting zones, the total mass m and volume V at any 
instant of time are given by 


m - m. + m 
b u 


(6) 


V 



(7) 


where p^and p u are spatially averaged densities of the burnt and unburnt gas 
mixtures. 

With reference to Figure 10, we also have 


m = m f + m^ = m f + p V;p (8) 

V = V, + Vq, (9) 

f f 

where the subscripts f and f denote the volumes inside and outside the flame 
front. The entrainment speed u g is defined as the average over the flame front 
surface of the normal component of the relative velocity at which unburnt gas 
mixture is entrained into the enflamed region 


p A,. u 
K u f e 


m f = - m f = Vf 


Cf 


( 10 ) 


Similarly, the average normal component of the unburnt gas speed u^ at the 
flame front is defined by 
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e u A f 


gf 


= ITb 


'A, 


p V; - p A 7U 
M u f F u pf p 


( 11 ) 


where u is the upward piston speed and the time derivative is taken for a 

r 

control volume- bounded- by the cylinder walls, the moving piston face of area 
Ap^ and a fixed surface just ahead of the flame' front. Adding equations 
(10) and (11) and making use of equation (5) which defines the average normal 
speed u^ of the flame front, we obtain the familiar relation 


= u , 


gf 


( 12 ) 


The evaluation of u g using equations (10) or (11) and (12) requires only the 
knowledge of V^, A^ and p u as functions of time. From the analysis discussed 
in sections 3.1 and 3.2, we have an estimate for such quantities. In connec- 
tion with the reliability of the estimate of u g it is useful to investigate 
the behavior of u g ^ at early times of flame propagation. Substituting equa- 
tion (6) into (7), taking the time derivative and recalling that V = -A u , 

P P 

we obtain 


p u 


P 


u 



(13) 


which substituted into (10) and (11) yields alternative expressions for u Q and 

u gf . At early times, when the flame has not yet hit the piston face, A p ^ = A p 

and it is reasonable to assume p. /p, ~ p /p and V-? ~ V to obtain 

b b K u *u f 



m. 


P A,. 
K u f 



( 14 ) 


Since u^ and u^ are determined independently, the. fractional error in 

and u can be estimated from the relations 
e 
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and 


/su ' 2 

gf • 

. K + iv 2 

vi 

1 v l 


0 o 1/2 

6u e 

( (6u f r + (6u qf ) Z ) 

u e 

u f “ u gf 


(15a) 


(15b) 


where the errors in and are assumed to be negligible. When the flame 
front first becomes visible, its radius is approximately 5 mm and the mass 
fraction burnt ~ 0.001 and its time derivative is essentially unmeasurable. 
Thus at early times the errors in u g and u^ are due entirely to the errors in 
x^ which, in practice, are large for x^ < 0.05. Nevertheless, a tentative 
estimate of u g at early times is obtained from Figure 9b where u^ and u^ are 
plotted for the typical cycle. The difference between the two curves yields 
the value of u . plotted in Figure 9c versus 6, corresponding to the integral 

c 

and direct estimates of u f . When x b > 0.05, such an estimate is reliable 

and in Figure 11 the mean values of u g are determined for two typical runs. 

When x^ <, 0.05, the reliability is weak as discussed above. The tendency of 

u to start from a low value (close to the value of the laminar flame speed S„), 
0 

then rapidly increase during the early period of combustion and eventually 

remain approximately constant for the remaining time is in agreement with the 

prediction of the rapid distortion theory (Wong and Hoult [29]). 

Table 5 shows the statistics of u g (average for x^ > 0.05), 

u^ = (R - x s )6N/A8 b , = A6 d /6N and the apparent ignition delay time • 

t. = (0. - 0 )/6N where 0. is defined in Figure 9a. No strong correlation is 
IIS i 

observed between u g and u^. Figure 12 shows the expected scaling of u g with 


engine speed. 
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3.4. Comparison of Enflamed and Burnt Mass Fractions and of Entrainment and 
Burning Rates 

The observation that for > 0.05 the entrainment speed u g does not show 
either increasing or decreasing trends is used to reconstruct flame front geo- 
metry in the non-visible region beyond the window edge. The dashed curve in 
Figure 9b is used to interpolate between the observed values of u^ and the 
limiting value of u^ = u g when the gas is fully entrained and u^ = 0. Knowing 
at any instant of time, the mass fraction enflamed x^ = m^/m can be esti- 
mated from equation (8) and normalized curves versus crank-angle for the 
sequence of Figure 4 are shown in Figure 13. It is observed that the mass 
fraction burnt curves are delayed with respect to the mass fraction enflamed 
curves of about 10 crank-angle degrees. Figure 14 shows normalized plots of 
x ^ vs. x^ for the same sequence. A significant amount of unburnt mass is ob- 
served inside the enflamed region (e.g., for the average cycle when x f = 0.6 
we read x b ~ 0.3). 

Figure 15a shows plots of normalized entrainment rate x f /x b max , burning 

rate V*b max and V*b max as a Unction of x b /x b max for the typical cycle. 
The importance of these plots is that they prove that the entrainment and 
burning processes are decoupled and thin flame assumptions cannot incorporate 
this experimental evidence. It is suggested that such plots be considered as 
important tests of turbulent flame propagation models. 

3.4.1. Characteristic burning time 

A characteristic burning time x can be defined as the instantaneous mean 
life time of an unburnt fluid particle in the enflamed region by the relation 

T x b - x f - x b . (16) 

To evaluate x from the available experimental information, either equation (16) 
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can be used or its time derivative 


T x* b = x f - x b (l + t ) (17) 

where the value for x can be estimated from the result of equation (16) (when 
x b > 0.2 we find x — 0.2). 

Figure 15b shows normalized plots of - x b and x^ - x b - By inspection 
of Figures 15a and 15b and using equations (16) and (17), both estimates of 
the value of x are shown in Figure 15c for the typical cycle, x is of order 
of 1 ms and after an early increase it shows a permanent decreasing trend. 

The notion of a characteristic burning time x was first introduced by 
Blizard and Keck [30]. It is to be noted that in the present work no use has 
been made of physical modeling interpretations of the analyzed phenomenon. 

3.4.2. Flame stretch factor 

At given thermodynamic conditions the values of the laminar flame speed 
are available from experimental measurements (e.g., Metghalchi and Keck [31]). 
Therefore, an equivalent laminar flame front area A b can be defined by the 
relation 

P„ A b s j - (18) 


Equation (18) substituted into (14) yields 






09) 


where the term A b /A-. is defined in the literature [32-35] 
as the flame stretch factor. Its estimate by using equation (18) is shown in 
Figure 9b. Also shown are the theoretical values at which the quantities 
shown would tend in the limit as r^ + 0 at 0 = 6^ if it is assumed A b /A^ 1 

and u £ •+■ S , i.e. , 
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1 

i 


(20a) 


u f " U 1 


= S 


£ p. 


(20b) 


where u-j and U 2 are the familiar relations of laminar flame propagation at 
constant pressure. The data would better agree with the assumption, say, of 
A fa /A f -v 1. 5-2-0 as r f 0. 


4. DISCUSSION AND CONCLUSIONS 
4.1. Previous Works 

As we repeatedly said, no physical modeling assumptions or interpretations 
have been used so far in the analysis of the experimental information. This 
is of fundamental importance for the results to be considered in testing physi- 
cal models and computer simulations of turbulent flame propagation in a spark- 
ignition engine. 

Blizard and Keck [30] have estimated u g using ionization probes in the 
place of photography. They obtained an average value of u g during a cycle, 
determined by assuming: (1) spherical flame geometry, centered at the spark 

plug; (2) an exponential distribution of burning times for the entrained 
eddies; and (3) constant u g along the cycle. A comparison of their results 
with the present estimates is given in Table 6. In their modeling work, 

Blizard and Keck defined a laminar burning length scale by the relation 

*1 - (21) 

interpreting it as the characteristic size of laminar subregions or "eddies" 
of the turbulent flow field which immediately start burning after entrainment 
at the laminar flame speed. A similar interpretation is given in flame 
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stretching theories [32-35] for the equivalent laminar 

flame front area as the total area of the extremely thin, highly convoluted 
laminar flame front contained within the turbulent burning zone behind the 
turbulent flame front. 

Lancaster, et al . [36] attempted to estimate u from the thermodynamic 

e 

information only. In fact, they estimated an equivalent thin-flame entrainment 
speed which is a good estimate only during the central period of combustion 
(around = 0). During the central period, using measured values of the 
turbulence intensity u 1 [37], they could correlate the values of u' , u g and 
for different operating conditions. 

The value of u g , in the literature more often called turbulent flame 
speed, has been related in a variety of ways to the turbulence intensity u' ■. 
and the laminar flame speed S^. This literature has been reviewed recently by 
Andrews, Bradley and Lwakabamba [38]. The most used relations are those found 
in [7-10]. A common feature of all of them is a direct dependence of u g 
on u’; i.e., if u 1 increases, u g also increases. The difference -between 
various models is in the nature of such dependence. 

Experimental measures of turbulence intensity and structure in the com- 
bustion chamber of an S.I. engine have also been presented in the literature 
[37, 39-43]. However, such measures, of fundamental impor- 

tance to understand the role of turbulence in the mechanisms of flame propaga- 
tion, are very critical and of difficult correct interpretation. A typical 
value of u' , for operating conditions similar to the present and arount TDC, 

is of the order of 10 m/s, i.e., of the same order as u . 

0 

4.2. Conclusions 

An experiment was carried out on a transparent-piston spark-ignition 
engine to obtain high-speed flame photography and cylinder pressure 
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measurements >of the turbulent flame propagation and combustion period of an 
engine operating cycle. 

1. - From the thermodynamic analysis of the pressure measurements only, 

data for cycl.e-to-cycle variations of both delay and burning angles 
are found to be of the same order of magnitude (about 20 percent). 

2. A new two-parameter empirical burning law., of practical usefulness 
in first-order computer simulations of engine operating cycles, has 
been presented and compared to experimental results. Its advantage 
i-s that it starts at spark and can simulate both delay and burning 
periods. Its application to simulate cyclic variations has been 
suggested. 

3. A model Independent method to analyze the combined photographic and 
pressure information is presented. Its importance is that the pre- 
sent results can be used to test physical modeling assumptions and 
■the results of computer solutions of the flow field equations. 

4,. The method is based on the definition of the 'entrainment speed u e ,, 
the uriburnt gas speed at the flame front u and the flame front 
speed u.p in terms of quantities that can be estimated from the experi- 
mental data. 

5. Estimated values of u g tend to increase rapidly in the early period of 
flame propagation until the mass fraction burnt is x^ « '0.05-0.10. 

Then u g remains approximately constant through the rest of the period. 
Its absolute value for the present engine operating conditions is 
about 10 m/s with cycle-to-cycle variations of the order of 20 percent. 

'6. Enflamed and burnt mass-fraction curves are compared through- the 

single cycle and a substantial fraction of the enflamed mass is found 
to be unburnt.. Thin-flame models cannot incorporate this observation. 
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7. Entrainment and burning rate curves are compared and the result 
proves that the two processes are decoupled. 

8. A characteristic burning time t is defined and found to assume values 
about 1 ms with a tendency to decrease during the combustion period. 

A flame stretch factor A^/A^ is defined and estimated from experimen- 
tal data. 
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APPENDIX: CHARACTERISTIC LENGTH SCALES 


To help the physical interpretation of the present results. Table A shows 
••a summary of characteristic length scales of the turbulent reacting flow field 
inside the combustion chamber. Columns a and b refer to the beginning and the 
.end of a typical cycle. Following the relations reviewed by Andrews, et al . 
,'[38‘] , we recall the expressions for Taylor and Kolmogorov microscales: 


where. 


•A = £ C 


1/2 Re/ 1 ' 2 


(A.1) 


n 


X 15 


-1/4 



(A. 2) 


l - integral length scale (assumed equal to the instantaneous 
clearance height) 

C - 48 (experimental constant for isotropic turbulence) 

U 1 J2. ' 

Re. = — — = turbulent Reynolds number 

Xj V 

U 1 A. 

Re. = — — - •= microscale ’Reynolds number 

A V 

u 1 = turbulence intensity (assumed - u & ~ 10 m/s) 
v = kinematic viscosity. 


£f = xu g can be Interpreted as a characteristic turbulent flame thickness. 

A diffusion length over the compression stroke is calculated from, the 
relation 6 = Alt , where l is the time from inlet valve closes. 

An 'estimate of'the laminar flame thickness, <S ~ k/(p S c ), where k 'is 

Xa U p 

the thermal conductivity, is also tabulated. 6^ is the value of the wall 
quench distance [23-24]. 
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TABLE 1.. Engine parameters 


Compression ratio- 

CR 

7.86 

Bore;, mm 

b 

101.6 

Stroke., mm 

S 

88'. 9- 

Clearance tie-ight, mm 

h 

13..0 

Equivalent wedge angle-,, rad 

a 

0.1.7 

Spark plug distance from center, mm 

x s 

-29. .5 

Inlet port diameter, mm' 

D 

48'.3 

Inlet- val*ve Tisft,. mm' 

L 

9. .3 

Inlet valve opens, deg-. ATDC 

I VO 

-k 

o 

o 

Inlet- valve cToses-, deg,. ATDC 

IVC 

I 

^ 1 

o 

o 

Exhaust valve opens,, de,g. ATDC 

EVO 

+T30° 

Exhaus.t valve closes-,, deg. ATDC 

EV.C 

+390° 

Spark- advance, deg. BTDC 

SA 

30° 

Throttle-, wide open 

WOT 


Fuel, iso^-octane 

C 8 H 18 


FueT-a.ir equivalence ratio 


1..0 - 1.2 

Crankshaft speed, rev/mi n 

N, 

900 -- 1200 



TABLE 2. Test matrix 


Run No. 

N 


1 

870 

1.09 

2 

1045 

1.00 

3 

1100 

1 .14 

4 

1120 

1.07 

5 

1220 

1.14 

6 

1230 

1.07 


N - Engine speed, rev/min 
$ = Fuel -air equivalence ratio 



TABLE 3. Statistics of A0 d , A0 b and A0 d + AG^ for the analyzed sequences 


Run 

No. 






Cycle No. 






Avg 

Rel . 
Var. 

First 

Cycle 



2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 


A 0 d 

20.0 

22.0 

19.1 

22.7 

20.4 

15.1 

18.5 

21.8 

19.6 

- 

- 

- 

19.9 

12% 

11.9 

1 

Ae b 

24.0 

19.9 

17.8 

22.6 

19.3 

25.8 

23.0 

28.3 

20.8 

- 

- 

- 

22.4 

14% 

20.3 


+ 

44.0 

41.9 

36.9 

45.3 

39.7 

40.9. 

41.5 

50.1 

40.4 

- 

- 

- 

42.3 

9% 

32.2 


A0 

25.2 

25.9 

33.8 

25.3 

32.5 

24.3 

28.8 

28.3 

22.1 

- 

- 

- 

27.4 

14% 

16.4 

2 

Ae b 

33.5 

24.1 

24.4 

29.3 

35.0 

35.5 

22.4 

23.4 

27.7 

- 

- 

- 

28.4 

17% 

23.2 


+ 

58.7 

50.0 

58.2 

54.6 

67.5 

59.8 

51.2 

50.7 

49.8 

- 

- 

- 

55.6 

11% 

39.6 


A 0 d 

20.9 

20.5 

19.1 

20.8 

19.0 

24.9 

28.1 

32.4 

20.0 

- 

- 

- 

22.9 

20% 

20.9 

3 

AG b 

18.2 

23.0 

29.9 

18.4 

21.9 

22.3 

27.8 

23.1 

24.1 

- 

- 

- 

22.6 

16% 

22.2 


+ 

39.1 

43.5 

49.0 

39.2 

40.9 

47.2 

55.9 

55.5 

44.1 

- 

- 

- 

46.0 

14% 

43.1 


A0 d 

24.8 

18.7 

20.7 

23.3 

23.0 

25.8 

16.8 

24.7 

22.2 

22.6 

17.5 

- 

21.8 

14% 

21 .2 

4 

A0 b 

18.3 

22.5 

22.6 

25.5 

22.1 

24.5 

30.4 

26.6 

19.6 

22.8 

29.0 

- 

24.0 

15% 

21.5 


+ 

43.1 

41.2 

43.3 

48.8 

45.1 

50.3 

47.2 

51 .3 

41 .8 

45.4 

46.5 

- 

45.8 

7% 

42.7 


A0 

22.6 

23.8 

30.9 

24.2 

19.1 

19.6 

22.1 

21.8 

24.6 

29.6 

24.9 

- 

23.9 

15% 

20.2 

5 

A6 b 

22.7 

20.4 

26.2 

19.6 

29.8 

24.8 

19.8 

44.3 

26.7 

28.7 

22.2 

- 

26.0 

27% 

19.6 


+ 

45.3 

44.2 

57.1 

43.8 

48.9 

44.4 

41.9 

66.1 

51.3 

58.3 

47.1 

- 

49.9 

15% 

39.8 


A0 

24.3 

30.2 

23.4 

20.8 

23.6 

23.5 

21.2 

23.8 

24.9 

25.4 

’24.8 

26.1 

24.3 

10% 

20.6 

6 

A9 b 

19.5 

31.6 

17.2 

30.4 

32.8 

28.9 

25.6 

20.4 

26.2 

25.1 

18.3 

23.8 

25.0 

21% 

18.8 


+ 

43.8 

61.8 

40.6 

51.2 

56.4 

52.4 

46.8 

44.2 

51.1 

50. '5 

43.1 

49.9 

49.3 

12% 

39.4 


TABLE 4. 


Summary of angles corresponding to data in Figure 4 


Cycle No. 

A0 d 

A0 b 

. A9 0.1 

A9 0.1~0.9 

A6 t 

"b" 

** 

1 

20.9 (0.94) 

22.2 

20.6 (0.93) 

22.5 (1.01) 

56. (2.52) 

0.40 

2 

20.9 (1.15) 

18.2 

19.8 (1.09) 

18.0 (0.99) 

50. (2.75) 

0.47 

3 

20.5 (0.89) 

23.0 

22.5 (0.98) 

21.8 (0.95) 

54. .(2.35) 

0.39 

4 

19.1 (0.64) 

29.9 

22.9 (0.77) 

28.8 (0.96) 

64. (2.14) 

0.30 

5 

-20.8 (1.13) 

18.4 

21.1 (1.15) 

18.3 (0.99) 

50. (2.72) 

0.46 

6 

19.0 (0.87) 

21.9 

21.0 (0.96) 

20.9 (0.95) 

54. (2.47) 

0.38 

7 

24.9 (1.12) 

22.3 

24.3 (1.09) 

22.4 (1.00) 

60. (2.69) 

0.46 

8 

28.1 (1.01) 

27.8 

23.9 (0.86) 

25.3.(0.91) 

62. (2.23) 

0.43 

9 

32.4 (1.40) 

23.1 

30.0 (1.30) 

23.4 (1.01) 

68. (2.94) 

0.55 

10 

20.0 (0.83) 

24.1 

21.2 (0.88) 

25.6 (1.06) 

62. (2.57) 

0.37 

Ave rage 

22.9 (1.00) 

22.6 

23.0 (1.01) 

22.7 (0.98) 

58.2 (2.54) 

0.42 


[0.99] 


[1.05] 

[1.02] 

[2.53] 


Rel . Var. 

20% (21%) 

16% 

13% (16%) 

15% (4%) 

11% (1055) 

17% 


* k 

Units: 0 in crank-angle degrees; u g in m/s. 

** 


First cycle not considered in averages. 




TABLE 5. 


Statistics of u g , = (R + r g JCN/AG^ in m/s 
and T d = A6/6N, x i = (6- - 0 $ )/6N in ms 


Run No. 







Cycle No. 






Avq 

Rel . 
Var., 

First 
C.ycl e 


2 

3 

4 

5 

6 

7 

8 

9 

10 

n 

12 

13 


Ug 

9.0 

10.6 

8.4 

9.1 


M 

— 

_ 

_ 


- 

- 

9.3 

10% 

8.4 


u b 

17.5 

21.1 

23.5 

18.5 

21.7 

16.2 

18.2 

14.8 

20.2 

- 

- 

- 

19.1 

15% 

20.6 

1 

T d 

3.83 

4.21 

3.66 

4.35 

3.91 

2.89 

3.54 

4.18 

3.75 

- 

- 

_ 

3.81 

12% 

2.28 


T i 

1.26 

0.68 

0.73 

0.62 

- 

- 

- 

- 

- 

- 

- 


0.82 

36% 

0.59 


u e 

8.4 

9.5 

8.0 

6.6 

9.0 

8.8 

11.8 

12.5 

11.9 

— 


— 

9.6 

21% 

9.5 

0 

u b 

15.0 

20.9 

20.6 

17.2 

14.4 

14.2 

22.5 

21.5 

18.2 

- 

- 

- 

18.3 

18% 

21.7 

u 

T d 

4.02 

4.13 

5.39 

4.04 

5.18 

3.88 

4.59 

4.51 

3.52 

- 

- 

- 

4.36 

14% 

2.62 


T i 

1.10 

0.87 

0.79 

0.52 

0.80 

1.32 

1.11 

1.40 

0.63 

- 

- 

- 

0.95 

32% 

1 .05 

> 


u e 

8.2 

8.3 

7.5 

8.4 

14,5 

10.6 

14.2 

13.0 

12.5 

— 

_ 

_ 

10.8 

26% 

8.4 

■3 

u b 

29.1 

23.0 

17.7 

28.8 

24.2 

23.8 

19.1 

22.9 

22.0 

- 

- 

- 

23.4 

15% 

23.9 

0 

T d 

3.17 

3.11 

2.89 

3.15 

2.88 

3.77 

4.26 

4.91 

3.03 

- 

- 

~ 

3.46 

20% 

3.17 


T i 

0.82 

0.81 

0.77 

0.86 

0.83 

0.76 

0.55 

1.05 

0.72 

- 

- 

- 

0.80 

17% 

0.66 


u e 

7.2 

14.8 

14.8 

13.0 

12.8 

9.6 

9.5 

12.1 

14.2 

13.8 

9.1 

_ 

11.9 

22% 

9.1 

A 

u b 

29.5 

24.0 

23.9 

21.2 

24.4 

22.0 

17.8 

20.3 

27.5 

23.7 

18.6 

- 

23.0 

15% 

25.1 


T d 

3.69 

2.78 

3.08 

3.47 

3.42 

3.84 

2.50 

3.68 

3.30 

3.36 

2.60 

- 

3.25 

14% 

3.15 


T i 

0.58 

0.37 

0.29 

0.32 

0.77 

0.52 

1 .03 

0.44 

0.00 

0.82 

0.71 

- 

0.53 

55% 

0.54 


u e 

9.7 

9.6 

10.8 

9.2 

9.9 

10.3 

15.6 

8.5 

12.2 


_ 

- 

10.6 

20% 

12.6 

Z 

u b 

25.9 

28.8 

22.4 

30.0 

19.7 

23.7 

29.7 

13.3 

22.0 

20.5 

26.5 

- 

23.9 

21% 

30.0 

0 

T d 

3.09 

3.25 

4.22 

3.31 

2.61 

2.68 

3.02 

2.98 

3.36 

4.04 

3.40 

- 

3.72 

15% 

2.76 


T i 

1.44 

0.68 

0.80 

0.47 

1.00 

0.46 

0.86 

0.62 

0.69 

- 

- 

- 

0.78 

39% 

0.58 


u e 

17.5 

9.4 

15.3 

11.7 

11.0 

11.9 

14.5 

9.3 

10.7 

11.7 

10.2 

15.5 

12.4 

21% 

12.3 

c 

u b 

30.4 

18.8 

34.5 

19.5 

18.1 

20.5 

23.1 

29.0 

22.6 

23.6 

32.4 

24.9 

24.8 

21% 

31 .5 

0 

T d 

3.29 

4.09 

3.17 

2.82 

3.20 

3.18 

2.87 

3.22 

3.37 

3.44 

3.36 

3.54 

3.30 

10% 

2.79 


T i 

0.81 

1.28 

0.46 

1.80 

0.78 

0.82 

0.97 

0.18 

0.91 

0.86 

0.47 

1.13 

0.87 

48% 

0.51 


TABLE 6. 


Entrainment speed, comparison with other parameters 


and Blizard and 

Keck’s results 




Present 

Blizard 



Work 

& Keck 

Engine speed, rev/mi n 

N 

1000 

1000 

Piston face inlet valve area ratio 

o. 

5.7 

8.4 

Mean piston speed, m/s 

S P 

3.0 

2.6 

Mean inlet gas speed, m/s 

^i 

17 

21 

Entrainment speed, m/s 

u e 

10 

6 

Rati o 

U e /D p 

3.3 

2.3 

Ratio 

u e /5 i 

0.59 

0.28 




TABLE A. Characteristic length sc lies 


Quantity 

Symbol 

Units 

"a" 

"b" 

pressure 

P 

MPa 

2 

4 

unburnt gas temperature 

T u 

K 

800 

920 

integral length scale. 

z 

mm 

13 

14 

turbulence intensity 

u 1 

m/s 

10 

10 

dynamic viscosity 

p 

mPa»s 

342 

378 

uriburnt gas density 

p u 

kg/m 3 

9.2 

16 

kinematic viscosity 

V 

? 

mm /s 

3.7 

2.4 

turbulent Reynolds number 

Re £ 


35000 

58000 

Taylor microscale 

X 

mm 

0.48 

0.40 

microscale Reynolds number 

Re x 


T300 

1670 

Kolmogorov microscale 

n 

pm 

6.8 

5.0 

Kolmogorov time scale 

T 

TJ 

ps 

12 

TO 

characteristic burning time 

X 

ms 

1.5 

0.9 

laminar flame speed 

S A 

m/s 

1.0 

1.5 

laminar burning- length scale 

l Z 

mm 

1.5 

1.3 

entrainment speed 

u e 

m/s 

TO' 

10 

flame "thickness 11 

£ f 

mm 

15 

9 

time after IVC 

l 

ms 

20 

23 

diffusion- length 

6 

mm 

0.28 

0.23 

laminar flame thickness 


pm 

5.4 

2.1 

wall quench distance 

6 q 

pm 

70 

40 
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•Figure 2. Pressure vs. crank-angle 0 for the typical sequence of cycles. 
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Figure 3. Typical photographic sequence corresponding to pressure curve 
number 9 of Figure 2. 
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Figure 4. Normal ized' mass fraction burnt n K x./x-. vs. e for pressure 
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curves in Figure 2. 





Figure 5. Definition of characteristic angles. 
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.Figure 6. f(n) = x b /x b max vs. n « x b /x b max . 





figure 7. Effect of parameters on the empirical burning law. 
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Figure 9a. Flame radius + x c - x g vs. crank-angle e for the typical contours in Fig. 

Figure 9b. Flame front speed evaluated with both direct and integral method and gas 
speed u gf vs. crank-angle for the same cycle. The flame stretch factor A fa /A f is 
also shown. 

Figure 9c. Entrainment speed u g vs. crank-angle. 
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Figure 10a. Cyl inder head geometry and assumed best spherical flame fronts for an early 
and a late visible contours. 

Figure 10b. Crank-angle or time vs. distance plot showing typical flame front and gas 

particle average paths. A geometrical visualization of u f , u^^ and. u g is also given. 
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Figure 11. Entrainment speed u g vs. for two typical sequences. 
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Figure 12. Entrainment speed u g vs. engine speed N. 
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Figure 13. Normalized mass fraction enflamed curves for the sequence of Figure- 2 to 4. 
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Figure 14. Enflamed vs. burnt mass fractions for the same sequence of Figure 13. 
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Figure 15a. Normalized rates of entrainment x,/x. and burning x./x K for the typical 
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cycle of Fig. 8. The second time derivative x*./x K has units of 1/ms. 
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Figure 15b. Normalized difference between enflamed and burnt mass fractions (x,.-x u )/x u 

f by b max 

in ms and normalized difference between entrainment and burning rates (x^-x. )/x. 
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Figure 15c. Characteristic burning time x vs. x,/x, calculated with two different 
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methods for the typical cycle. 



